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ABSTRACT 

We use high resolution images obtained with the Advanced Camera for Surveys on board the 
Hubble Space Telescope to determine morphology, nuclear luminosity and structural parameters of the 
spheroidal component for a sample of 20 Seyfert galaxies at z = 0.36. We combine these measurements 
with spectroscopic information from the Keck Telescope (paper I) to determine the black hole mass 
- spheroid luminosity relation (Mbh-Lr), the Fundamental Plane (FP) of the host galaxies and the 
black hole mass - spheroid velocity dispersion relation (Mbu-cf). The FP is consistent with that of 
inactive spheroids at comparable redshifts. Assuming pure luminosity evolution, we find that the 
host spheroids had smaller luminosity and stellar velocity dispersion than today for a fixed Mbh- 
The offsets correspond to AlogLB,o = 0.40 ± 0.11 ± 0.15 (A log Mbh = 0.51 ± 0.14 ± 0.19) and 
A log cr = 0.13 ± 0.03 ± 0.05 (A log Mbh = 0.54 ± 0.12 ± 0.21), respectively for the Mbh-L and Mbh-ct 
relations (the double error bars indicate random and systematic uncertainties). A detailed analysis 
of known systematic errors and selection effects shows that they cannot account for the observed 
offset. We conclude that the data are inconsistent with pure luminosity evolution and the existence 
of universal and tight scaling relations. In order to obey the three local scaling relations by z = - 
assuming no significant black hole growth - the distant spheroids have to grow their stellar mass by 
approximately 60% (AlogMgph = 0.20 ± 0.14) in the next 4 billion years, while preserving their size 
and holding their stellar mass to light ratio approximately constant. The measured evolution can be 
expressed as Mnn/Msph cx (1 + z)^-^^^-^, consistent with black holes of a few 10^ Mq completing their 
growth before their host galaxies. Based on the disturbed morphologies of a fraction of the sample 
(6/20) we suggest collisional mergers with disk-dominated systems as the physical mechanism driving 
the evolution. 

Subject headings: black hole physics: accretion — galaxies: active — galaxies: evolution — quasars: 



general 



1. INTRODUCTION 

In the local Universe, most g:alactic nuclei harbor a 
supermassive black hole (e.g.. iKor mend v fc Richstond 
[1995). The mass of the black hole correlates with global 
properties of the host, such as the velocity dispersion 
and luminosity of the spheroidal (or bulge) component, 
indicating a connection betwee n nuclear activity and 
galaxy fo rmation and evolution (Magorrian et al. 1998"; 



Ferrarese Merritt 2000; Gebhardt et al. 2000, 2001; 

Hunti>2 003: HaringjLRi2dL2004; Novak et al. 
Driver 2Q07.). Understanding the ori- 
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gin of this relation is a major challenge for cosmological 
models and is believed to hold the key to solving several 
astrophysical problems such as the role of feedback from 
nuclear activity in suppressing star formation in massive 
galaxies (e.g.. Granato et al. 2004: Di Matteo et al, 200j.; 



Croton et al 1 120061 : iMalbon et al ]l2006l : ICiotti fc Ostriked 



2003). 



In the standard cosmological scenario, spheroids grow 
by mergers of smaller galaxies while black holes grow 
by accreting surrounding matter. Depending on the rel- 
ative timing of the two processes, the scaling relations 
between black hole mass and spheroid luminosity and 
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velocity dispersion (hereafter Mbh-L - or Mbh-Lb to 
specify the blue band - and Mbh-c" relations, respec- 
tively) could also evolve with cosmic time. For example 
- if the spheroid evolved passively due to aging of stellar 
evolutions, changing L but not a, while the black hole 
grows by accretion - we would expect evolution in both 
the M-L and the M-<j relations, with the latter evolving 
more slowly than the former. By contrast, if black holes 
completed their growth first and the dominant mode of 
growth now is the transformation of stellar disks into 
spheroi ds, the two relations would evolve in the opposite 
sense (Croton 2006). The tightness of the local rela- 
tionships has been interpreted as evidence for feedback 
that synchronizes the relative growth. In this context, 
the tight relationships would be naturally reproduced 
if galaxies and black holes moved up the M-cr relation 
during merging events, so that the correlation would ap- 
pear not to evolve with redshift (Haehnelt fc Kauff mannI 
M). 

Detailed theoretical predictions are extremely difficult 
due to the daunting range of scales - ranging from the 
Mpc scale halo to the pc scale dynamical sphere of influ- 
ence of the black hole, to the /ipc scale of the accre- 
tion disk - and physical processes involved- radiative 
transfer, heating and cooling, accretion, just to name 
a few. In spite of the challenge, numerous groups have 
been able to develop models that are increasingly more 
successful at reproducing a variety of observation s (e.g . 
[Granato et al. 2004; Miral da-Escud e fc Kollmei ej l2005l: 
iHopkins et al., ,2006b : .Croton et all 12006. : .Malbon et all 
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[2QQ6I : iDe Lucia fc BlaizQtl[2QQ7l : iCiotti fc Ost7ik^[2QQ7h . 
However, to this date, the evolution of scahng laws re- 
mains very uncertain and sensitive to the schemes and 
approximations adopted to deal with the c omplex physics 
(coni pare for example the recent works bv lHopkins et alJ 
I2QQ7I and [Robertson et al. 2006 ). 

Accurate empirical measurements are needed to dis- 
criminate between scenarios, and provide input on the 
relative importance of various physical phenomena as 
well as on the accuracy of approximations. With this 
goal in mind, a number of groups have started ob- 
servational programs to trace the evolutio n of scaling 
laws o yer cosmic time (IShield s et al.l 120031; 'Treu et al. 
200^ IWalteretal.1 [200^; [Ad elberger & Steidel 2005; 
Woo et al.'2006: McL ure et al.ll2006i :iPeng et al. 2006a|lb|; 
Salviander et al. 2006^. However, observers face two 
fundamental limitations. On the one hand, - since 
the sphere of influence of supermassive black holes is 
too small to be resolved at cosmological distances with 
present technology - black hole mass estimates can only 
be obtained for active galaxies. Typically, this involves 
broad line AGN and the dynamics of the broad line re- 
gion, with consequent loss of information about the host 
galaxy which is swamped by nuclear light. On the other 
hand - at a deeper level - the evolution of the scaling 
laws, depends on the interplay of at least four physical 
processes: i) black hole accretion; ii) evolution of the stel- 
lar populations; iii) dynamical evolution of the spheroid, 
e.g. through mergers; iv) black hole feedback on star 
formation. For this reason, even when a scaling law can 
be measured as a function of redshift, the interpretation 
is often times amb ig;uous. For example, as discussed by 
iPeng et all (|2006bl ) , how much of the evolution of the M- 
L relation is due to evolution in the spheroid mass and 
how much is due to evolution of the stellar populations? 

In this paper we address the two fundamental limita- 
tions by adopting the following strategy. First, as in our 
pilot study (.Treu et aT]l200j ) and in the first paper of 
this series (|Woo et al.l I2006L hereafter paper I), we fo- 
cus on relatively moderate redshift {z ~ 0.36) and lumi- 
nosity (monochromatic luminosity at 5 100 A ~ 10^^ erg 
s~^) AGN. Although the lookback time is considerably 
smaller than that of the most distant quasars studied 
by other groups, this choice allows us to determine the 
host galaxy properties with considerably smaller uncer- 
tainties. Second, we concentrate on a relatively small 
sample and measure several independent properties of 
the host galaxies. Paper I reported stellar velocity dis- 
persion measurements based on Keck spectroscopy. This 
paper presents host spheroid luminosity and size mea- 
surements based on HST-ACS imaging. We use this com- 
bined dataset to study at the same time and for the same 
sample the Mbh-L and M-<j relations and the Fundamen- 
tal Plane of host spheroids. The combination of these 
diagnostics - which can be thought as projections of a 
more fundamental manifold (e.g., Hopkins et al. 2007) - 
allows us to disentangle stellar mass growth, stellar popu- 
lation evolution and black hole gro wth, and helps to iden- 
tify the processes at work (see also ICiotti fc van Albadal 
120011 : iNipoti et al.ll2003l : iBovlan-Kolchin et al.l[2006h . 

The paper is organized as follows. In § O we summa- 
rize the properties of the sample, describe the HST-ACS 
observations, and present our surface photometry. We 



also derive black hole masses using nuclear luminosities 
as measured from HST images and the new calibration 
of th e broad line region size - nuclear luminosity rela- 
tion (|Bentz et al.l [2006b). together with H/3 line widths 
from paper I. To construct a suitable local comparison 
sample, we use Sloan Digital Sky Survey images to derive 
new measurements of the spheroid luminosity of a sample 
of local S eyferts with Mrh ava ilable from reverberation 
mapping ([Peterson et aP 120041 : iBentz et al.ll2006a[ ). §[3] 
describes our main results, i.e. the Fundamental Plane, 
Mbh-Lb, and Mbh-c" relation of the distant Seyferts. 
Detailed estimates of systematic errors and selection ef- 
fects are given in § [H § [5] analyzes the three scaling 
relations under the assumption that the distant Seyferts 
will evolve to match the local relations and derives con- 
straints on the evolution of stellar mass, size and stellar 
populations of the host spheroids as a function of black 
hole growth. The results are discussed and compared 
with the literature in § [6l and summarized in § [71 

Throughout this paper magnitudes are given in the AB 
scale. We assume a concordance cosmology with matter 
and dark energy density Qrn = 0.3, = 0.7, and Hubble 
constant Ho=70 kms~^Mpc~^. 

2. DATA 

2.1. Sample selection and observations 

The sample discussed in this paper coincides with that 
discussed in paper I of this series, with minor excep- 
tions. Two of the objects (S16 and S31) studied in this 
paper do not have Keck spectroscopy due to the unfa- 
vorable weather conditions at Keck when follow-up ob- 
servations were planned. One of the objects with deep 
Keck spectroscopy (S28) has no ACS images because it 
was observed after the schedule of the Hubble program 
had been completed. As in our previous paper, except 
for S99 that was select ed before SDSS data release 1 
(DRl) (|Treu et al.ll2004f ). the objects were selected from 
the SDSS archive, based on the following criteria: i) spa- 
tially resolved in the SDSS images; ii) redshift between 
0.35 and 0.37; iii) H/3 equivalent width and Gaussian 
width greater than 5 A in the rest frame. Most of the 
objects were selected from DR-1, and the remaining ones 
(S21,S23, S26, S27) were selected from DR-2. After the 
initial selection based on these criteria, the SDSS spec- 
tra of the objects observable from Keck were visually 
inspected by two of us (TT and MAM) to check line iden- 
tification. Objects showing strong Fe II nuclear emission 
(the main obstacle to velocity dispersion measurement) 
were eliminated from the sample. For example, for the 
run of September 2003, out of 33 observable objects, 8 
were rejected on the basis of visual inspection (the mean 
and rms r' magnitude of the selected and rejected sam- 
ples are consistent, 18.76 and 0.42 vs. 18.98 and 0.63, 
respectively). Of the remaining 25 objects, 12 were ob- 
served during the run, based on observability. The frac- 
tion of rejected is similar for following runs, although the 
total number of objects increases significantly in subse- 
quent data releases. No color selection was imposed, al- 
though a post-facto analysis shows that the SDSS colors 
of the sample are intermediate between those of a quasar 
and those of an old stellar population, consistent with 
the comparable fraction of nuclear and stellar light in- 
ferred from the HST images in the rest of this paper. 
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Fig. 1. — Postage stamp images of the 20 Seyfert galaxies in the sample. Each postage stamp is 12 arcseconds on a side, i.e. approximately 
60kpc. Merging (Interacting) objects are identified by the label M (I). The centers of galaxies Sll, S31 and S99 are obscured by dust lanes, 
which prevented accurate surface photometry of the nucleus and spheroid. 

the Fundamental Plane. This section presents the mea- 
surements with an extensive discussion of systematic and 
random uncertainties. 



Coordinates, redshifts and other basic properties of the 
sample are given in Table [H 

The sample was observed using the Wide Field Camera 
of the ACS on board HST between August and December 
of 2004 as part of General Observer program 10216 (PI: 
Treu). Each object was observed for one orbit, split in 
four exposures dithered by semi-integer pixel shifts to en- 
sure cosmic ray and defect removal as well as to improve 
sampling of the point spread function (PSF). The total 
exposure times range between 2148s and 2360s. Filter 
F775W {i') was chosen so as to sample the region corre- 
sponding to the rest frame 5100 A in order to estimate 
the size of the broad line region, and to avoid contamina- 
tion from the broad emission lines H/3 and Ha. This filter 
choice also provides spheroid luminosity redward of the 
4000 A break, and it is thus appropriate to infer spheroid 
structural parameters and rest frame B and V luminosity 
with minimal uncertainty due to filter transformations. 
Postage stamp images of the targets are shown in Fig- 
ure [TJ 

2.2. Reduction and analysis 

The following observables are needed to investigate our 
science questions: i) spheroid luminosity; ii) nuclear lu- 
minosity, to estimate the size of broad line region, and 
hence the black hole mass; iii) effective radius and ef- 
fective surface brightness of the spheroid, to construct 



2.2.1. Reduction 

The images were reduced using multidrizzle to re- 
move cosmic rays and defects, correct for distortion and 
improve sampling of the PSF. Based on Monte Carlo 
simulations we adopted a final pixel size of 0V04, which 
provides the best compromise between sampling of the 
PSF, signal to noise ratio of the individual pixels, and 
noise correlation. 

2.2.2. Surface photometry 

Surface photometry was derived using the GALFIT soft- 
ware (Peng et al. 2002) to fit two-dimensional models to 
the data. Optimizing a large number of non-linear pa- 
rameters is a notoriously difficult problem, so we pro- 
ceeded in steps, adding one component at a time. First, 
we fitted a point source to determine the center of the 
system, which is assumed common to all components. 
We th en added a spheroid modeled as a ide VaucouleursI 
('1948') r"*^/^ profile, and/or an exponential disk, if re- 
quired by the statisti c. If needed, an additional com- 
ponent described by a Sersic' (1968) profile with index 
n=0.5 was added to model the bar, identified as an elon- 
gated feature with a change in position angle with respect 
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Fig. 2. — Surface photometry I. The surface brightness profile measured from the data are shown together with that determined from 
the two-dimensional model that best fits the data. The surface brightness profile of each component is shown separately to illustrate the 
relative contribution as a fun ction of radius. Note that this plot is for illustration only, and the fits were performed in twodimensions 
as described in Section 12.2.21 using cutouts of 20^^ on a side. The early truncation of a few profiles (e.g. S09) is only an artifact of the 
elliptical isophotes routine used to make the plots, due to nearby objects. For the measurement, the nearby objects were modeled and 
fitted simultaneously in twodimensions. 



to the overall surface brightness distribution (see, e.g., 
SOI in Figure 1). Bars were found to be present in 7/20 
cases and contribute between 4% and 16% of the total 
light. Extensive tests were conducted to ensure that the 
solution corresponded to the true global minimum of the 
over the parameter space. Bright unsaturated stars 
in the field of the images were use to create a library of 
43 point spread functions. Each galaxy was fitted with 
all the PSFs in the library to find the best fitting one as 
well as the best fitting parameters. 

In three cases (Sll, S31, S99), a prominent dust lane 
prevented us from determining accurate surface photom- 
etry. In two cases (803 and SIO) no stable solution with a 
sizable spheroid could be found, as the spheroidal compo- 
nent tended to become vanishingly small in size. There- 
fore we fixed the spheroid half light radius to 2.5 pixels 
(0.1 arcsec), the minimum size that could be resolved 



given our PSF and we considered the measured lumi- 
nosity as an upper limit. Similarly, for three objects 
(S12, S21, S23) the measured spheroid half light radius 
is very close to our resolution limit (i.e. <3pixels), and 
therefore we also consider their luminosity as an upper 
limits. Extensive testing shows that the upper limit to 
the spheroid luminosity is robust with respect to the 
choice of the fixed half light radius, and of the PSF. 
Twelve out of twenty cases provided stable best fitting 
models and hence robust measurements. For illustration 
purposes one-dimensional surface brightness profiles (ob- 
tained with IRAF task ellipse) are shown in Figures [2] 
and [31 

The formal statistical uncertainties on the spheroid lu- 
minosity are typically 0.05 mags. The rms scatter of 
the parameters for all the statistically acceptable PSFs 
was adopted as the systematic uncertainty due to PSF 
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Fig. 3. — Surface photometry II. As in Figure [2] for the rest of the sample. 



modeling. We estimated the total uncertainty - includ- 
ing systematic errors - by varying systematically all the 
fitting parameters and PSF, finding that the results are 
typically stable within 0.2 mags. Spheroid luminosity 
was then transformed into rest frame B-band luminosity 
by correcting for Galactic extinction and a pplying K- 
color c orrections, calculated as described in iTreu et al.l 
(|2QQlbi ). Errors on extinction and K-color corrections 
are a few hundredths of a magnitude. Conservatively, 
we adopt 0.5 mags (i.e. 0.2 dex) as the total uncertainty 
on spheroid luminosity. This uncertainty is smaller than 
the estimated uncertainty on black hole mass from sin- 
gle epoch measurements ~0.4 dex and therefore sufficient 
to meet our goal of constructing the M-Lb relationship. 
Observed and intrinsic spheroid luminosities are given in 
Table El 

Another source of systematic uncertainty is the choice 
of the surface brightness profile used to describe the 
spheroid. Although r^/"^ profiles are the traditional and 
widely used choice in the analysis of AGN host galax- 



ies (e.g. iBentz et al.l l2006al ). detailed studies of nearby 
bulges show that Sersic profiles with lower Sersic index 
(typically 2-3) can provide a better fit. To estimate sys- 
tematic errors on the spheroid luminosity associated with 
our choice of profile, we repeated our analysis using a Ser- 
sic profile with index 2 and 3 (Figured)). As expected, we 
find that the best fit spheroid luminosity decreases with 
the Sersic index. Quantitatively, the spheroid luminosity 
decreases on average by 0.15 (0.33) mags when changing 
the Sersic index from 4 to 3 (2). Similarly, the point 
source luminosity increases on average by 0.02 (0.05) 
magnitudes and the effective radius of the spheroids de- 
creases by 9% (11%) when changing the Sersic index from 
4 to 3 (2). We thus conclude that adopting a r^/^ profile 
provides a conservative measurement of the maximum 
luminosity of the spheroid and of the minimum nuclear 
luminosity (and hence Mbh)- We will come back to this 
point in the discussion of our results in Section 14. 2[ 

Other measurements of interest are total magnitude, 
nuclear luminosity at rest frame 5 100 A (L5100), the 
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Fig. 4. — Systematic effects due to adopted spheroid profile. 
Changing the Sersic index of the spheroid from n=4 to n=3 (2), 
reduces systematically the spheroid luminosity by 0.15 (0.33) mag- 
nitudes, and increases the nuclear luminosity by 0.02 (0.05) mag- 
nitudes. 

combination of effective radius (Re) and effective sur- 
face brightness (SBe) that enters the Fundamental Plane 
FPph = log Re - 0.32 SBe (see [31]) and the fraction of 
nuclear light (fnuc)- The total uncertainties, conserva- 
tively estimated as for the spheroid luminosity, are 20% 
on L5100 and on fnuc, and 0.1 on FPph. The relevant 
quantities are also listed in Tabled 

2.3. Black hole mass 

As in paper I of this series, black hole masses 
were derived from the width of H/3 and the ob- 
served nuclear luminosity at 5100 A, using the so- 
called 'virial' metho d or empirical l y cali b rated photo- 
ioniz a tion method (IWandel et al] 1 19991 : IVestergaardI 
I2OO2I : [V^stergaard & Peterson" "2006). Briefly summa- 
rized, the method assumes that the kinematics of the 
broad line region trace the gravitational field of the cen- 
tral black hole. The width of the line provides the ve- 
locity scale, while the nuclear continuum luminosity pro- 
vides the size via the empirical correlation observe d for 
the local reverberation mapp ed samp le (Wandel et al.l 
I1999I : iKaspi et al.ll200d l2005l : iBentz et al.ll2006b[ ). The 
virial coefficient is obtained by requiring that local AGN 
hosts and qui escent galaxies obey the same M-a relation 
(|Qnken et al.l l2Q04l : |Greene & Hd l2006[ ) . This method 
has been shown to provide an estimate of the black hole 
mass within a factor of 2-3. For this paper we will as- 
sume as in the previous paper of the series a nominal 
uncertainty of 0.4 dex on black hole mass obtained with 
this method (see Peterson 2007, for a recent discussion 
of the method). However, since we allow for unknown in- 
trinsic scatter in the scaling relations when fitting for the 
intercept (§[3]), adopting a larger random uncertainty on 
each individual black hole m ass estimate (e.g. 0.6 dex, 
IVestergaard fc Petersonll2006l ) has a negligible impact in 
terms of overall uncertainty (§ 14. 3p . 

As in paper I, we adopt the second moment of the 
broad component of H/3 as our velocity scale since this 
is more robustly measured than the alternative FWHM 



(paper I; see also Peterson et al. 2004; Collin et al. 2006). 
The velocity scale is measured from Keck spectra when 
available and from SDSS spectra in the few cases when 
Keck spectra are not available (see paper I for details). 
As shown in paper I, our methodology gives unbiased 
black hole masses when applied to the local sample of cal- 
ibrators. A major improvement with respect to paper I 
is that we can use high resolution Hubble images to mea- 
sure the nuclear luminosity and infer broad line region 
size. In practice we use the following equation to estimate 
black hole masses, obtained combining the most recent 
calibratio ns of the size luminosity relation and v irial co- 
efficients (|Onken et al.ll2004l : iBentz et al.ll2006aD . 



logMBH = 8.58+2 log 



V3000kms" 



+0.518 log 



5100 



lO^^ergs- 
(1) 

where a^p is the second moment of the broad H/3 line 
profile. We note that for the sample in common with 
paper I, the use of nuclear luminosities and this new cal- 
ibration implies black hole masses smaller by 0.09 dex on 
average, consistent with our estimate of the magnitude 
of this systematic error in paper I. Black hole masses are 
listed in Table 

2.4. Local comparison samples 

To study evolutionary effects we consider two compar- 
ison samples to define the local M-Lb relation. First 
we consider the s a mple of quiescent galaxies collected by 
iMarconi fc HuntI (|2003l ). The sample consists of galax- 
ies with black hole mass determined from spatially re- 
solved kinematics and with spheroid luminosity deter- 
mined from two-dimensional surface photometry. The 
second sample is that of local Seyfert galaxies with black 
hole mass measured via reverberation mapping discussed 
in Peterson et al. (2004) and Onken et al. (2004). To en- 
sure self-consistent determination of spheroid luminosity 
for the local and distant Seyfert sample we measured the 
parameters of the local sample in exactly the same way as 
we do for distant galaxies. We searched the SDSS archive 
for g'-band images of local Seyferts, which provide a very 
good match to our distant galaxies in terms of resolution 
(the seeing is typically 10 times the HST PSF, but the an- 
gular size distance is typically 10 times smaller) and rest 
frame wavelength. SDSS images for nine Seyferts were 
found in the archive. The resulting spheroid luminosi- 
ties for the local Seyfert samples are listed in Table [3l 
together with redshifts, black hole masses and velocity 
dispersions from the literature. 

The two samples are complementary in terms of vices 
and virtues. The first sample is larger in size and black 
hole masses, and spheroid luminosities are most robustly 
measured: i) Mbh is obtained from spatially resolved 
kinematics; ii) the determination of spheroid luminosity 
does not suffer from poor resolution or the presence of 
a prominent point source in the center. However, the 
second sample is the most appropriate for a direct com- 
parison for a variety of reasons: i) the virial coefficient 
is unknown, but assuming that it is not evolving with 
redshift, the comparison between Seyfert samples is in- 
dependent of its numerical value; ii) if the local M-L rela- 
tion is not universal, Seyferts may define a different M-L 
relation than quiescent early-type galaxies - perhaps be- 
cause they are at a different evolutionary stage. Thus, 
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we conclude that at this stage the more conservative ap- 
proach is to consider both samples for the following anal- 
ysis and consider the uncertainty on the local relation as 
an additional source of systematic errors. More compre- 
hensive studies of the local M-L relation are needed to 
eliminate this source of uncertainty. This is left for future 
work. 

3. RESULTS 

This section presents the main results of this paper. 
First in § 13.11 we discuss the Fundamental Plane of the 
host galaxies in comparison with that of normal quiescent 
galaxies, for the subsample of objects that have both a 
and structural parameters. In § 13 21 we present the Mbh- 
Lb relation for distant Seyfert galaxies. Having estab- 
lished that the Fundamental Plane of the host galaxies 
is indistinguishable from that of quiescent galaxies, we 
adopt the luminos ity evolution inferred from FP analysis 
(|Treu et al.ll2QQ5al jbl) to compare with the local Mbh-Lb 
relation. In § 13.31 we revisit the Mbh-cj" relation derived 
in paper I, using the new improved black hole mass es- 
timates based on nuclear luminosities determined from 
HST imaging. 

3.1. Fundamental Plane 

In the local Universe early-type galaxies obey the Fun- 
damental Plane, i.e. an empirical correlation between 
effective radius, central velocity dispersion and effective 
surface brightness of the form: 

logi?e = c.loga + /3SBe+7, (2) 

with a = 1.25, (3 = 0.32 and 7 = -9.00 for the Coma 
Cluster in the B(AB) band. The evolution of the FP 
out to redshift z ~ 0.4 is well established for quies- 
cent early-type galaxies (e.g., iTreu et aI]|2005U l2001a . 
and references therein) , for AGN hosts (Woo et al.""2004", 
I2OO5), and for the spheroidal component of spiral galax- 
ies with bulge-to-total luminosity ratio greater than 0.2 
(MacArthur et al. 2007, in preparation). 

In Figure [5] we plot the FP parameters of the sub- 
set of distant Seyferts for which structural parameters 
and stellar velocity dispersion is available, together with 
a comparable sam ple of quiescent galaxies taken from 
iTreu et al] (|2005al lbl). The good agreement between the 
FP of quiescent and active galaxies gives us confidence 
that our surface photometry is not systematically biased 
by the presence of a nuclear point source. The offset with 
respect to the local relation (solid line) is normally in- 
terpreted as evolution of the stellar populations. Under 
the assumption of pure luminosity evolution, luminos- 
ity evolves with redshift such that the expected value at 
z = is given by 

log Lb,o = log Lb - (0.72 ± 0.06 ± OM)z (3) 

(jTreu et al.ll2005b[ ). Applying the same assumption to 
the sample of distant Seyferts at z = 0.36 implies that: 
log Lb, = log Lb — 0.26 ± 0.03, where random and sys- 
tematic errors have been added in quadrature for sim- 
plicity. 

3.2. Mbh-Lb Relation 

In order to compare with the local relations we need 
to account for the fact that the luminosity of a stellar 
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Fig. 5. — Edge on view of the Fundamental Plane of the dis- 
tant Seyfert galaxies (circles) and of normal quiescent galaxies at 
comparable redshift (squares). Seyfert galaxies for which only an 
upper limit to the spheroid luminosity is available are plotted as 
cyan symbols. The errors on the FP variables are highly correlated, 
and points are allowed to move mostly within the plane. The error 
component perpendicular to the plane is estimated to be ~ 0.1 dex 
for the Seyfert sample when projected along the log Re axis. Most 
of the apparent thickness of the FP for quiescent galaxies is due 
to the relatively large range in redshifts and luminosity evolution 
during the corresponding interval in cosmic time. The local FP of 
the Coma Cluster is shown as a solid black line for comparison. 

population decreases as it ages. As our first working 
hypothesis, in this Section we will present our results as- 
suming that stellar populations evolved as inferred from 
the FP studies under a pure evolution scenario, and use 
the variable Lb,o obtained using equation [3l A more 
general discussion will be given in § [5l 

The relation between black hole mass and host 
spheroid luminosity for our points as well as for local 
comparison samples is shown in Figure [H The main re- 
sult is that the Seyfert galaxies at z = 0.36 and those at 
z = cover approximately the same range in spheroid 
luminosity, but the average black hole mass is higher 
for the distant Seyferts. The mismatch is exacerbated 
if one considers that five of the distant Seyferts mea- 
surements are upper limits (while only one of the lo- 
cal Seyferts is an upper limit). Conservatively we will 
generally consider the measured offset as the best esti- 
mate of the offset although it should be kept in mind 
that it is most likely a lower limit. Quantitatively, 
the offset with respect to the fiducial local relationship 
(solid line) corresponds to ALb,o = 0.32 ± 0.11 ± 0.15 
(A log Mbh = 0.40 ± 0.14 ± 0.i9), considering the in- 
trinsic scatter of the relation to be a free parameter and 
then marginalizing over it. Listed systematic errors are 
derived as discussed in Section [H If measured with re- 
spect to the local Seyferts (magenta line), the best es- 
timate of the offset changes to A log Mbh = 0.63. This 
is visualized in Figure [71 where we plot the distribution 
of residuals with respect to the local fiducial relation for 
the distant and local samples. 

We conclude that pure luminosity evolution is inconsis- 
tent with the observations, if we required that at z = all 
galaxies obey the Mbh-Lb relation. We cannot solve this 
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Fig. 6. — Black hole mass spheroid luminosity relation. Distant 
Seyferts are shown as large circles with error bars. Blue circles 
represent measurements, cyan circles with leftward arrows repre- 
sent upper limits to the spheroid luminosity. Local Seyferts are 
shown as large squares with error bars. Black hole masses are from 
iPeterson et al. (2004) and Bentz et al. (2006b), spheroid lummosi- 
ties are from this work (Table [3]). Magenta squares represent mea- 
surements, the red square represents an upper limit to the spheroid 
luminosity. Small black squares represent local quiescent galaxies 
from (Marconi & Hunt 2003). The black line is the best fits to 
the Mar coni & Hunt (2003) sample. The magenta line is the best 
fit to the local Seyferts. For direct comparison the spheroid lu- 
minosity of the distant Seyfer ts h as been evolved to z = using 
logLB,o = log Lb - 0.26 (see §[32]). 

inconsistency by invoking different luminosity evolution. 
The stellar populations would be required to become 
brighter with time in order to reconcile the data with 
the local relationship. This can be ruled out on physi- 
cal grounds. Therefore, taking this result at face value, 
we have to conclude that either not all spheroids obey 
the M-L relations, or that a significant amount of new 
stars have been added to the spheroids since z = 0.36. 
The interpretation of this result will be discussed in Sec- 
tions [5] and [6l after we conclude presenting the evidence 
and discussing systematic errors and selection effects in 
the remainder of this section. 

3.3. Mbh-ct Relation 

Figure [8] shows the Mbh-c" relation for the distant 
Seyferts as well as the local comparison samples. The 
samples and symbols are the same as in Figure [6] with 
few minor exceptions: i) only the 14 distant Seyferts 
with available stellar velocity dispersion (from paper I) 
are plotted, including S28 and S99 for which HST pho- 
tometry is not available, ii) local relation s are from 
iTremaine etHI (j2002|) and iFerrarese fc FordI ([2005); iii) 
Local Seyferts obey the same relation as quiescent galax- 
ies by construction, as this is t he constrai nt used to derive 
the virial coefficient ( Onk en et al.ll2063 ), so there is no 
need to show a separate local relation for Seyferts. The 
only substantial difference with respect to Figure 8 in 
paper I is that black hole masses have been recalculated 
based on HST photometry. 

The Mbh-<^ relation shows the same basic result as 
the Mbh-Lb relation. At fixed host galaxy properties. 



Fig. 7. — Distribution of residuals in logMsH with respect to the 
fiducial local relation of Marconi & Hunt (2003, solid black line in 
Figure [6]). The Upper panel shows distant Seyferts (measurements 
in blue, upper limits in cyan; note that upper limits in logLB,o 
correspond to lower limits in A log Mbh)- The lower panel shows 
the distribution of residuals for the local samples of quiescent and 
Seyfert galaxies. 




(7(km s~^) 

Fig. 8. — Black hole mass velocity dispersion relation. The large 
blue circles represent distant Seyferts. Black hole masses are from 
this work, velocity dispersions are from paper I. The large magenta 
squares represent local Seyferts, the small black squares represent 
local quiescent galaxies. The solid and dashed lines are the best fit 
relations from Tremaine et al. ( 2002) and Ferrarese & Ford ( 2003). 
The local Seyferts obey the same Mbh-<^ relation as local quiescent 
galaxies by construction, see Onken et al. (2004) for discussion. 

distant Seyferts appear to host larger black hole mass 
than local ones. Quantitatively, the offset corresponds 
to Alogcr = 0.13 ± 0.03 ± 0.05 (i.e. A log Mbh = 
0.54±0.12±0.21). Note that offset is slightly reduced 
(0.08 in A log Mbh) with respect to paper I and the 
systematic error is smaller due to the improved esti- 
mate of black hole masses that avoid stellar contami- 
nation. As shown in Figure [8] an d discussed in paper I , 
adopting the local relation from lTremain"e et all (|2002r ) 
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or iFerrarese fc FordI (|2QQ5l ) changes the offset by much 
less than the error bars because the two local relations 
are very well matched in the range of black hole mass 
and stellar velocity dispersion covered by our sample. 

4. SYSTEMATICS AND SELECTION EFFECTS 

This section is devoted to understanding and estimat- 
ing systematic errors. First, in § 14.11 and § 14. 2[ we list 
potential sources of systematic error and estimate as 
accurately as possible the associated uncertainty. For 
brevity we will not repeat the analysis of systematics 
already discussed in paper I, unless when substantial 
changes/improvements are introduced. § 14.31 derives the 
uncertainty on our estimate of the random errors due to 
our assumed accuracy of single-epoch black hole mass 
determinations. Then, in § 14. 4[ we introduce a Monte 
Carlo scheme to simulate unknown selection effects and 
estimate potential biases. § 14.51 gives a very short sum- 
mary of this Section for the impatient reader. 

4.1. Is Mbh overestimated? 

Understanding the Mbh estimate is clearly important. 
While velocity dispersion and spheroid photometry are 
measured independently and present mutually consistent 
results (as supported by the FP relation, which does not 
involve Mbh), only one measurement of Mbh is available 
and that determines the evolution of both relationships. 
Our measurement of Mbh depends on two observables 
(velocity scale and nuclear luminosity), on the empiri- 
cally calibrated relationship between nuclear luminosity 
and broad line region size and on the choice of the virial 
coefficient. We now go through each source of error and 
estimate the systematic error that they may introduce 
for a sample of our size. 

4.1.1. Hp line width and velocity scale 

In paper I we demonstrated that our measurement of 
the velocity scale from the second moment of H/3 in sin- 
gle epoch spectra is unbiased. This means that - when 
applied to the sample of local calibrators - our estimator 
yields the same Mbh as obtained by the original studies 
based on the second moment of 11/3 m easured fr o m the 
variable part of the spectra. Recently, iPetersonl (|2QQ7l) 
suggested that single epoch spectra may overestimate the 
line width by ~20%, corresponding to an error of ~0.15 
dex on black hole mass, if the same virial coefficient is 
assumed. This effect is likely to depend on the setup of 
the experiment and on the variability pattern. However, 
to be conservative, we will consider this as a global un- 
certainty to the zero point of our measured black hole 
mass. I Woo et al.l (j2QQ7l ) studied the time variability of 
the width of H/3 for a subset of our objects. They found 
that the r.m.s. variations of the second moment of 11/3 
are less than 14% (i.e. ~ 0.1 dex in Mbh)- For our 
sample of 14-17 objects this has therefore negligible ef- 
fect on the average with respect to the aforementioned 
source of error. Our best estimate of overall systematic 
uncertainty due to line width measurement is 0.15 dex. 

To check for systematic effects due to systematic 
change of the line shape with redshif t or sampl e prop- 
erties such as Eddington Ratio (e.g.. iCollin et al. 1120061 ) 
we also computed black hole masses using the full width 
half maximum instead of the second moment of the line 



as velocity scale. After removing the narrow compo- 
ne nt of H/3 u s ing; [O III] as a temp l ate a s described 
in iTreu et all (|2004 f) and IWoo et all (|2006i ). we mea- 
sured the FWHM (McGih et al. 2007, in preparation) 
and derived black hole mass estimates Using Eq. 1 in 
iNetzer fc Tr akhtenbrot (2007). This alternative esti- 
mate of black hole mass (Mbh,nt) agrees very well with 
that adopted in this paper, with an average offset of 
(logMBH/MBH,NT) = -0.02 ±0.05 dex, i.e. consistent 
with no offset. 

4.1.2. Nuclear luminosity and size luminosity relation 

The new calibra tion of the size-luminosity relation 
(|Bentz et al.ll2006a[ ). and our HST based nuclear lumi- 
nosities are a substantial improvement with respect to 
paper I, allowing us to avoid issues related to host galaxy 
contamination. The dominant source of error is now 
the intrinsic scatter of the size luminosity relation (30%) 
corresponding to ~ 0.1 dex in Mbh per individual ob- 
ject. For our sample size, this translates into a negligible 
source of error in the mean. As far as the overall calibra- 
tion of the relation is concerned, for the luminosity range 
of our sample, the 68% range on the slope and intercept 
of the relationship translate into an overall uncertainty in 
size of order 0.02 dex. We thus conclude that the uncer- 
tainty due to the (known) scatter of the size luminosity 
relation is negligible. 

4.1.3. Virial coefficient and zero point of the local relation 

The final source of error is related to the virial co- 
efficient needed to convert size and velocity scales into 
mass and thus to the zero point of the local relation. As 
discussed in detail in paper I, we adopt the virial coef- 
ficient determined bv'Onke n et al.l (|2004 !) requiring that 
the M-(T relation be the same for local quiescent and ac- 
tive galaxies. The uncertainty on the average of the virial 
coefficient is 34%, thus corresponding to an uncertainty 
of 0.13 dex on the calibration of the black hole mass. 
However - as discussed in paper I - assuming that the 
virial coefficient does not evolve with redshift, this fac- 
tor cancels out between the local and distant samples 
and therefore is irrelevant. Incidentally, we note that 
our results can also be interpreted as cosmic evolution 
of the virial coefficient, if one is willing to drop this as- 
sumption. We consider this an unlikely explanation, as 
it would require the geometry or kinematics of AGN to 
evolve with cosmic time, but unfortunately it will not be 
possible to discard it until direct measurements of the 
virial coefficient can be obtained in some other way. 

The situation is slightly more complex for the Mbh- 
L relation. Once the virial coefficient is fixed there is 
no more freedom. Therefore the difference between the 
intercept of the relation for local active and quiescent 
(0.23 dex in Mbh) is an additional source of uncertainty. 
Larger samples of galaxies with well determined black 
hole mass and spheroid luminosity are needed to over- 
come this limitation. For the purpose of this paper we 
consider the difference as an additional source of uncer- 
tainty. To produce a single measurement, we weight the 
two samples equally (for reasons discussed in § 3) and 
take the average and semi-difference of the two intercepts 
as best estimate of the local average and systematic un- 
certainty on the zero point (~ 0.12 dex). 
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4.2. Is Lb underestimated? 

The main sources of uncertainty that could affect the 
spheroid luminosity of the sample as a whole are system- 
atic errors in K-corrections, in the adopted luminosity 
evolution, and in the choice of the surface brightness pro- 
file for the bulge. The uncertainty on the K-correction is 
negligible, at most 0.02 dex, as estimated using a range 
of stellar population models to compute the transforma- 
tion from observed F775W to rest frame B. The luminos- 
ity evolution as ni easured from the Fundamental Plane 
(jTreu et al.ll2005aD carries an overall uncertainty of 0.03 
dex. Even assuming that this translates to a shift of the 
whole sample, this is still a negligible source of error. As 
discussed in Section 12.2.21 adopting a Sersic profile in- 
stead of a r^/^ profile tends to systematically reduce the 
bulge luminosity while increasing the point source lumi- 
nosity, and therefore black hole mass, thus moving points 
further away from the local M-L relation of quiescent 
galaxies. In contrast, since we used the same identical 
technique for the local Seyferts, adopting a Sersic profile 
would not alter the observed offset. We conclude that 
the measured offset is not overestimated by an apprecia- 
ble amount due to known systematic errors on spheroid 
luminosity determination. 

4.3. Are random errors underestimated? 

Finally, the uncertainty of black hole mass estimates 
from sing;le epoch data is assumed to be 0.4 dex 
(|Vestergaard^ l2002l ) . This represents all sources of ran- 
dom error that contribute to the scatter in Eq. [TJ After 
quantifying sources of systematic errors on the mean of 
EqlU we now consider whether the random error associ- 
ated with our measurement of the offset of the scaling re- 
lations could be underestimated. To test this, we repeat 
the analysis assuming that the r andom error on each indi- 
vidual Mbh estimate is 0.6 dex ("Vestergaard & Peterson 
[2OO6) . The inferred random error on A log Mbh increases 
by 0.03 dex, e.g. from 0.14 dex to 0.17 dex for the Mrh- 
L relation. Note that our analysis includes a nuisance 
parameter to account for the unknown intrinsic scatter 
of scaling relations and therefore the estimated errors are 
slightly different than what would be naively derived con- 
sidering only errors on the y axis (0.4 / — 1 = 0.4/4 
= 0.1 dex vs 0.6/\/N — 1=0.15). This is a negligible er- 
ror on the error, considering that the systematic term is 
dominant. The total error, defined as the quadratic sum 
of the random and systematic uncertainties, would only 
change by 0.01 dex. 

4.4. Selection bias 

In this section we estimate possible bias due to selec- 
tion effects. Our galaxies are selected based on their nu- 
clear properties (having a broad line AGN). Therefore, 
intrinsic scatter and observational errors in the scaling 
relations could conceivably lead us to favor large black 
hole masses and therefore overestimate evolution^. 

Our present sample covers approximately a decade in 
black hole mass (10^ — lO^M©). The upper limit is natu- 
rally expected because of the steep drop of the black hole 
mass function. Is the lower limit in mass a result of some 

^ This bias is similar to Malmquist bias for luminosity selected 
samples of standard candles 



unknown selection effect? Our nominal selection limits 
on line width and flux are small enough that we would 
have been able to select objects with black hole masses as 
low as 10^ M0, as verified by running our measurement 
tools on the entire SDSS-DR5 spectroscopic database at 
this redshift. However the objects with Mbh well below 
IO^Mq are a small fraction of the total. This maybe be 
due to an intrinsic drop in the black hole mass function, 
to a selection effect in the SDSS spectroscopic sample or 
to a combination of the two. The decline in the number 
of objects at 2; = 0.36 with Mrh well below 10 ^ is also 
seen in Figure 1 of lNetzer fc TrakhtenbrotI (|2007l ). where 
most of the points lie above the group identified as hav- 
ing Mbh between 10^"^ and 10^*^ M©. Estimating the 
black hole mass function is beyond the scope of this work 
(e.g-. lBernardi et al. 2006; Greene fc Ho 2007). However, 
to understand the implications of this selection effect on 
our measurement, we can take a conservative approach 
and model this drop as a sharp selection in Mbh- In the 
following discussion we assume for reference that our im- 
plicit selection function log Mbh /M© > 7.9 (the results 
are unchanged if the limit is chosen to be 7.8 or 8.0). 
We then use Monte Carlo simulations to determine the 
amount of bias on the offset of the intercept of the scal- 
ing relations introduced by the implicit selection process. 
The scaling relations are populated according to the ve- 
locity dispersion fu nction and spheroid lu minosity func- 
tion d etermined bv ISheth et al.l (|2003f ) and iDriver et al.l 
(|2007l ) respectively, and errors on both axis are taken into 
account. We assume that the intrinsic scatter of the rela- 
tion is smaller than the measurement errors (0.4-0.5 dex) 
which is consistent with local estimates. Note that in this 
scheme all the selection procedure is modeled a single 
hard threshold in black hole mass and we only work in 
the scaling relation variables. This simple scheme allows 
us to bypass all the assumptions that would be needed to 
simulate from first principles the observational selection 
effects which are a complex combination of total flux, 
broad line flux, broad line width (at both ends), bulge 
luminosity and stellar velocity dispersion. 

The results of the simulations are shown in Figure [H 
The curves show the recovered ("measured") offset as a 
function of the input offset A log Mbh- The upper panel 
shows the results for the Mbh-<^ relation, the lower panel 
shows the results for the Mbh-L relation. The bias is 
almost negligible for the Mbh-c" relation, while it is at 
most 0.1 dex for the Mbh-L relation. The difference be- 
tween the two relations is due to: i) the smaller errors 
on velocity dispersion than those on spheroid luminosity; 
ii) the different behavior of the two distributions at the 
faint (low velocity) end. The velocity dispersion func- 
tion peaks in the range covered by our sample, while the 
luminosity function of bulges extends to much smaller lu- 
minosities. As a softer cutoff would imply smaller bias, 
we conclude that selection effects are responsible for at 
most 0.1 dex of the observed offset of the Mbh-L relation, 
and that this source of bias is negligible for the Mbh-c^ 
relation. 

4.5. Summary of uncertainties and best measurements 

In conclusion, the total systematic error on the evolu- 
tion of Mbu-o- relation is 0.21 dex in Mbh (i-e. slightly 
reduced with respect to paper I due to the improved nu- 
clear luminosity estimates), and 0.19 on the Mbh-L re- 
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Fig. 9. — Estimate of possible bias due to selection effects. The 
effects on an hypothetical selection log ( Mb h /M q) > 7.9 are cal- 
culated via Monte Carlo simulations (see ^ 14.41 for details). The 
measured offset is shown as a function of the simulated input off- 
set in Mbh with respect to the local relation. Upper panel: 
results for the Mbh-<7 relation. The bias is negligible due to 
the small error on a (compared to that on Mbh) and the flatness of 
the velocity dispersion function in the range of interest. The mea- 
sured offset and the corresponding input offset are identified by the 
solid blue lines. Dashed lines identify error bars obtained by sum- 
ming in quadrature random and systematic errors. Lower panel: 
results for the Mbh-Lb relation. The bias is somewhat more 
significant due to the larger error on Lb and to the steeper faint 
end of the spheroid luminosity function. The solid blue line and 
large filled circle identify the lower limit to the offset as measured 
with respect to the average of local quiescent and active galaxies. 
Dashed lines represent error bars as in the upper panel. 

lation, dominated by the uncertainty on the black hole 
mass and on the uncertainty on the local relation. 

Although we caution the reader to keep in mind all 
the caveats discussed above, we now condense all the 
information discussed in this section in seven numbers. 
Our best estimates of the offset of the relations - without 
accounting for selection effects are: 

• Mbh-L: a log Mbh = 0.51 ± 0.14 ± 0.19 

• Mbh-ct: a log Mbh = 0.54 ± 0.12 ± 0.21 

Unknown selection effects could remove as much as 0.1 
dex to the offset of the Mbh-L relation. 

5. EVOLUTION OF THE FUNDAMENTAL MANIFOLD OF 
BLACK HOLES AND SPHEROIDS 

In the previous sections we concluded that pure lumi- 
nosity evolution is inconsistent with the observations and 
the requirement that the distant Seyferts lie on the local 
relations. In this section, we explore the evolutionary 
constraints that can be obtained by examining the evo- 
lution of black holes and spheroids in a four dimensional 
parameter space - with axes given by luminosity, veloc- 
ity dispersion, size and black hole mass - instead that on 
lower dimension spaces. 

In other words, we know that in the local universe, 
black holes and their host spheroids lie on the FP, the 
Mbh-c" and Mbh-L relations - 3 relationships involving 4 
parameters. Thus, by requiring that our distant galaxies 
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Fig. 10. — Inferred evolution of the fundamental manifold param- 
eters (solid blue lines), as a function of change in black hole mass. 
Dotted blue lines represent the 68% confidence bands. The symbol 
6 indicates the difference between the parameter at z = 0.36 and 
the parameter at z = 0. The red dashed horizontal lines repre- 
sent pure luminosity evolution of the host galaxies for reference. 
If there is no black hole growth (5 log Mbh =0) the host galaxies 
of the distant Seyferts need to increase their velocity dispersion 
by 0.13 ± 0.06 dex and increase their luminosity by 0.20 ± 0.19 
dex in the next 4 Gyrs in order to obey simultaneously the local 
FP, M.^ii-a, and Mbh-Lr relations, with no significant change in 
the effective radius. For comparison, pure luminosity evolution of 
the spheroid would predict 0, -0.26 and 0, respectively (dashed red 
lines). This corresponds to a growth of ~60% of the spheroid stel- 
lar mass and approximately constant stellar mass-to-light ratio. If 
there is significant black hole growth (51ogMBH< 0) the increase 
in cr and Lb and the decrease in effective radius must be more 
pronounced. 

land on the local relationships at z = 0, we can derive 
constraints on the evolution in size, luminosity and mass 
of the host galaxies, as a function of black hole growth 
in the same time span. 

The results of this exercise are shown in Figure [TOl For 
simplicity, systematic and random errors are combined 
in quadrature in order to compute confidence bands, and 
neglecting covariance. The main result is that in order to 
satisfy all scaling relations in the local universe the host 
galaxies have to grow in luminosity and velocity disper- 
sion, while leaving the size substantially unchanged. 

Consider the no black hole growth scenario (consistent 
with the low Eddington ratios, see paper I). Assuming 
that Mgph cx (J^Re^ the spheroid mass has to increase by 
0.20 ± 0.14 dex (- 60%), by today. In other words. 



Ml 



BH 



Msph 



oc (1 + z) 



1.5±1.0 



(4) 



Similarly, the spheroid luminosity has to increase by 
0.20 ± 0.19 dex, implying that its average stellar mass 
to light ratio has to stay approximately constant (with 
an uncertainty of approximately 50%) over the next 4 
Gyrs. This requires an injection of younger stars to 
counteract the ageing of the resident population. The 
effective radius does not change significantly ^logi^e = 
0.07±0.07. Those trends are amplified if significant black 
hole growth is assumed. 

Although at this stage the measurement uncertainties 
do not allow firm conclusions, it is clear that this em- 
pirical methodology holds great promise as a way to 
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disentangle evolution of the stellar populations, black 
hole growth, and dynamical evolution of the host. In 
the future, with larger samples covering a wider range 
of masses and redshifts, it will be possible to study in 
the detail the co-evolution of spheroids and black holes, 
perhaps in the framework or more funda mental underly- 
i ng co rrelations such as that proposed by iHopkins et alJ 
(l2QQ7h . 

6. DISCUSSION AND CONCLUSIONS 

In previous sections, we have presented evidence that 
just four billion years ago, there was a population of su- 
permassive black holes living in smaller and less luminous 
spheroids than today. Based on a detailed analysis of 
uncertainties and selection effects, this evolutionary sig- 
nature appears significant at the 95% CL. In this section 
we will briefiy discuss our results in the broader context 
of galaxy formation and evolution and identify a possi- 
ble physical mechanism. However, before discussing the 
interpretation of this result, it is worth mentioning two 
caveats (see also the discussion in paper I), i) Our mea- 
surement relies on Mbh estimates based on an empiri- 
cally calibrated method. The sample of local calibrators 
is very small and does not necessarily match exactly the 
properties of our distant sample. More work on the local 
and distant samples remains to be done before unknown 
systematics can be firmly ruled out. ii) The local scaling 
relations are based on very limited samples. The qui- 
escent sample consists of approximately 40 objects. In 
addition to the intrinsic difficulty of resolving the sphere 
of infiuence, the local sample is composed of earlier type 
galaxies than the sample considered here. As discussed 
in paper I, the two samples may have an evolutionary se- 
quence built in at the selection. In this respect, the local 
AGN sample is very important, as the selection process 
and properties are very similar to those of our own sam- 
ple. However the local AGN sample is comparable or 
smaller in size to our distant sample, and therefore it is 
hard to identify subtle differences which could point to 
some unknown selection effect. 

Keeping these two caveats in the back of our minds, we 
move on to discuss and interpret the main result of this 
paper, that black hole mass growth appears to be com- 
pleted before bulge growth. An important clue is given 
by the morphology of the host galaxies, provided by the 
HST images. The majority of galaxies in our sample 
are not spheroid-dominated elliptical or lenticulars, but 
rather intermediate or late type spirals. Most remark- 
ably, as shown in Figure, 6/20 galaxies are observed to 
show signs of a major ongoing merging (3) or to be mor- 
phologically disturbed (3). Considering only the three 
most extreme cases and assuming that major mergers are 
visually identifiable for ~0.5-l Gyr (Co x et al. .2006h «a). 
this finding implies that most or possibly all the galaxies 
in our sample will undergo a major merger in between 
the time of observation and today. Gas-rich mergers are 
believed to be the main mechanism transforming stel- 
lar and gaseous disks into stellar spheroids, leading to 
substantial increase in spheroid luminosity and velocity 
dispersion. By converting rotation supported stars into 
pressure supported stars, a merger with a disk domi- 
nated system (and hence negligible supermassive black 
hole) could grow the bulge more efficiently then the cor- 
responding growth of the black hole by accretion of cold 



gas (jCrotonll2006h . 

As calculated in § O in order to satisfy the local scal- 
ing relations, assuming that black hole growth is neg- 
ligible, the stellar mass of the spheroid will have to 
grow by approximately 60%, while the spheroid mass 
to light ratio wold have to remain approximately con- 
stant, qualitatively cons istent with a single dissipative 
merger (|Cox et al.ll2006bH al). and associated bursts of star 
formation. This scenario is similar to that discussed in 
paper I, to which the reader is referred for further dis- 
cussion and references, with two important additions: i) 
the high resolution HST images provide direct evidence 
for gas rich mergers; ii) our joint analysis of several scal- 
ing relations allows us to quantify the expected growth 
of the spheroid mass and constrain the evolution of the 
stellar populations of the spheroid via the evolution of 
the mass-to-light ratio. Our imaging results may also 
suggest an apparent evolution of Seyfert 1 host galaxy 
morphologies from z = 0.36 to the current epoch. This 
is because most local Seyfert 1 galaxies have substantial 
bulges, with normal stellar po pulations, and no signs of 
interactions (jHunt et al.ll 19991 ) . Although currently most 
Seyfert Is at z=0 do not show signs of ongoing mergers, 
they do have hints tha t these occurred ~ 1 Gyr earlier 
(jHunt fc Malkanlll999l ). but that their host galaxy mor- 
pholoffles have s ettled back down to normal by today 
(jHunt Malk^[200l . So our result appears to be con- 
sistent with the id ea that Seyfer t host galaxies, as quies - 
cent galaxies (e.g. iBundv et al"]l2004l : iLotz et all [20061 ). 
were more disturbed by interactions 4 Gyrs ago than 
they are today. Unfortunately our sample is too small to 
make a proper comparison, and to derive statistically sig- 
nificant results. However, we checked the consistency of 
this statement using the GOODS database. The merger 
rate for our sample is 3/20 (O.I3I0 07) considering only 
close pairs/mergers, and it increases to 6/20 considering 
all disturbed systems (0.30 ± 0.12). As a control sample, 
we selected from GOODS all galaxies with luminosity 
17.6 < i' < 20.0 (i.e. the same range in stellar lumi- 
nosity that our sample, excluding the point source), in 
the redshift range 0.26 < z < 0.46 (morphologies, pho- 
tometric and spectroscopic redshifts from Bundy et al. 
2006) and performed the same visual classification. We 
found 8/42 close pairs/mergers (0.19 ± 0.07), increasing 
to 12/42 (0.28±0.08) if we consider ah disturbed systems. 
This is in good agreement with the fraction observed for 
the distant Seyferts, and somewhat larger th an observed 
in the local universe (e.g. lPatton et al.ll2002l ). The small 
size of our sample does not warrant a more detailed at- 
tempt to measure the merger fraction, so this is left for 
future work and larger samples. 

Independent evidence appears to support the scenario 
discussed above. From the galaxy evolution point of 
view, although the most massive spheroids appear to 
have completed most of their growth hy z ~ 1, it is 
clear that spheroids of a few 10^^ M© are undergoing 
significant evolution at 2: ~ 0.4. This view is sup- 
ported by the signatures of recent star for mation de- 
tected by Fundamental Plane studies (e.g., PlVeu et al.l 
2Q05bl : Ivan der Wei et al.l[2005l : [di Serego Alighieri et~all 
2005L and references therein), by the evolution of the 
mass function (e.g. Bundy, Treu & Elhs 2007 and refer- 
ences therein), and by the evolution of the quenching or 
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transition mass (e.g. jHopkins et aI]|2QQ7l , and references 
therein). From the point of view of black hole demo- 
graphics, it is hard to pinpoint with sufficient precision 
the growth of black holes in this mass range, although 
general arguments based on the global star formation 
and accretion history suggest that black hole growth may 
predate bulge growth (e.g., iMerlonil |2QQ4[ ). At larger 
masses, t he existence of impressively luminous quasars 
at z - 6 (iFan et all [20061 ) - with the high black hole to 
host galaxy mass ratio determined by radio observations 
(Walter et al. 2004), appears to be consistent with our 
scenario, although it is hard to make a direct compari- 
son, considering that evolution may very well be mass- 
dependent. 

Recent studies of the scaling relations between black 
hole mass and host galaxy properties tend to focus on 
higher m asses and redshifts than are in our sample. For 
example, iPeng et al] (|2006bl jal) study the host galaxies of 
lensed quasars out to 2; ^ 4, ruling out pure luminosity 
evolution and finding that the ratio between Mbh and 
Mgph was ~ 4 times larger at 2; ~ 2 — 3 than today. 
Similar results are obtained by othe r studies of quasar 
host galaxies at coni parable redshifts ([Shields et al.l2006l : 
iMcLure et aI]l2006D . These results are consistent with 
our own, although the comparison requires caution con- 
sidering that the higher redshift studies typically rely 
on UV broad emission lines and fiuxes for estimating 
Mbh, instead of and that the contrast between 
AGN and host light is more unfavorable than in our 
case. This prevents accurate decomposition of the host 
spheroid light and direct determination of the stellar ve- 
locity dispersion, which needs to be inferred from proxies 
such as CO and J O III] lin e widths (iBonning et al]l2005l : 
ISalviander et al.l 12006: W alter et al.ll2004|). An alt ernate 
approach followed by Adelberg er fc St eidell (|2005l ) leads 
to the opposite conclusion. They use the correlation 
length of AGN hosts at z ~ 2 — 3 to estimate the virial 
mass of the halo, and the C IV line width and UV fiux 
at 1350 A to estimate Mbh- They compare the inferred 
relation between halo mass and black hole mass with the 
local relation (Ferrarese 2002), finding no evidence for 
evolution. Given their error bars, they rule out evolution 
of over one order of magnitude in the ratio - i.e. evolu- 
tion of the form MBu/Msph cx (1 + z)^-^, with z ^ 2.5 - 
at 90% CL. However, they cannot rule out evolution by 
a factor of 6 that would be predicted extrapolating our 
best estimate. 

In conclusion, it seems that several lines of evidence 
are beginning to point in the same direction: black holes 
appear to complete their growth before their host galax- 
ies. The uncertainties are still large and much work 
remains to be done before evolution can be considered 
conclusively detected. However, with the advent of new 
technologies such as laser guide star adaptive optics and 
robotic telescopes it will hopefully be possible to improve 
dramatically our ability to measure black hole masses in 
the local and distant universe, and thus reduce the main 
source of error. 

7. SUMMARY 

This paper is devoted to the study of the cosmic evo- 
lution of black holes and their host galaxies. To this aim 
we have performed a detailed analysis of a sample of 22 
Seyfert 1 galaxies at z ~ 0.36. The choice of this par- 



ticular redshift and moderate luminosity AGNs allows 
us to derive precision measurements of the host galaxy 
properties, as well as to obtain an estimate of the black 
hole mass from the dynamics of the broad line region. 
Using high resolution images taken with ACS we derived 
luminosity, effective radius and effective surface bright- 
ness of the host spheroid as well as nuclear luminosity. 
We combined this information with emission line widths 
and stellar velocity dispersion based on high signal to 
noise Keck spectroscopy (paper I) to construct the Mbh- 
L, Mbh-c" and Fundamental Plane relations of distant 
Seyferts. We compared the z ~ 0.36 scaling relations 
with those followed by local samples of quiescent and ac- 
tive galaxies to determine evolutionary trends. The main 
results can be summarized as follows: 

1. The Mbh-Lb relation at z ~ 0.36 is inconsistent 
with the local relation and the assumption of pure 
luminosity evolution of the host galaxy. Adopt- 
ing pure luminosity evolution consistent with Fun- 
damental Plane studies, the offset from the local 
relation corresponds to an offset in present day B- 
band luminosity of AlogLB,o = 0.40 ± 0.11 ± 0.15, 
i.e. A log Mbh = 0.51 ± 0.14 ± 0.19, in the sense 
that black holes lived in smaller bulges at z ~ 0.36 
than today. 

2. The Mbh-c" relation at z ~ 0.36 is inconsistent with 
the local relation and the assumption of pure lumi- 
nosity evolution. The offset with respect to the lo- 
cal relation corresponds to A log a = 0.13 ± 0.03 ± 
0.05, i.e. A log Mbh = 0.54 ± 0.12 ± 0.21, in the 
sense that black holes lived in smaller bulges at 
z ~ 0.36 than today. 

3. Monte Carlo simulations show that the offset is not 
the result of selection effects, which are negligible 
for the Mbh-c" relation, and can account for at most 
0.1 dex of the observed offset of the Mbh-L relation. 

4. In order to satisfy the local Mbh-c^, Mbh-L and 
FP relations by z = - assuming no black hole 
growth - our distant spheroids have to grow their 
stellar mass by approximately 60% (AlogMgp = 
0.20 ± 0.14) in the next 4 billion years, while pre- 
serving their size and holding their stellar mass 
to light ratio approximately constant. This corre- 
sponds to an evolution of the black hole to spheroid 
mass ratio of the form Mbh/Msph cx (1 + z)^*^^^*^. 

Assuming that our results are not due to unknown 
systematic errors or unknown selection effects, the ob- 
served evolution can be qualitatively explained if our 
Seyferts undergo a single collisional merger with a disk- 
dominated system between z=0.36 and today. This is 
consistent with the observed merging/interacting frac- 
tion and a timescale for merging visibility of ~1 Gyr. A 
single merger could increase the spheroid mass by trans- 
porting stellar mass from the progenitors disks, without 
the corresponding growth of the central black holes due 
to the lack of black hole in the disk dominated system. At 
the same time, this process would add younger stars to 
the spheroid (either from the merging disks or from newly 
formed stars) thus counteracting the fading of the old 
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stellar populations and producing an approximately con- 
stant stellar mass to light ratio in the spheroid. Numer- 
ical simulations including realistic prescriptions for star 
formation, AGN activity and mass loss will be needed to 
see if these mergers do, indeed, preserve Re and Msph/Le- 
If these indications are supported by future studies, then 
they will confirm that black holes completed their growth 
before their host galaxies and are perhaps to be seen less 
as a by-product of galax y formation than as an orches- 
trator (jSilk Reeslll998l : iBlandfordlllQQOT ). 
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TABLE 1 
Sample properties 
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S28 


16 


11 


56.30 


+45 


16 11.04 


0.3682 


18.59 


210+10 


S29 


21 


58 


41.93 


-01 


15 00.33 


0.3575^ 


18.77 




S31 


10 


15 


27.26 


+62 


59 11.51 


0.3504^ 


18.14 




S99 


16 


00 


02.80 


+41 


30 27.00 


0.3690 


18.33 


224+12 



Note. — Col. (1): Target ID. Col. (2): RA. Col. (3): DEC. 
Col. (4): Redshift from stellar absorption lines. Col. (5): Extinction 
corrected i' AB magnitude from SDSS photometry. Col. (6): Stellar 
velocity dispersion in kms"-*^ from paper I. 
^ redshift from SDSS DR4. 



TABLE 2 

New measurements 



Name 


i' (total) 


i' (spheroid) 


log Lb/Lq^b 


SBe,B 


Re 


L5100 


fnuc 


logMBH/Mo 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


SOI 


18.50 


19.92 


10.28 


21.85 


5.29 


0.74 


0.29 


8.21 


S02 


19.03 


19.85 


10.31 


20.27 


2.63 


0.36 


0.22 


7.99 


S03^ 


17.94 


20.23 


10.16 


17.04 


0.50 


1.69 


0.39 


8.29 


S04 


18.06 


20.12 


10.20 


18.36 


0.96 


1.42 


0.36 


8.45 


S05 


17.93 


20.45 


10.07 


18.84 


1.03 


2.04 


0.47 


8.77 


S06 


18.35 


20.48 


10.06 


18.81 


1.01 


0.54 


0.18 


8.17 


SOT 


17.79 


20.35 


10.11 


18.69 


1.01 


2.26 


0.45 


8.55 


SOS 


18.31 


21.75 


9.55 


20.50 


1.23 


1.25 


0.40 


8.10 


S09 


18.17 


19.00 


10.65 


19.87 


3.24 


0.78 


0.22 


8.15 


SIO^ 


18.01 


19.30 


10.53 


16.08 


0.49 


1.11 


0.27 


8.27 


S12^ 


18.12 


21.16 


9.78 


18.14 


0.54 


1.05 


0.28 


8.69 


S16 


19.14 


22.26 


9.34 


19.97 


0.76 


0.70 


0.48 


8.26 


S21^ 


17.45 


18.95 


10.67 


15.82 


0.51 


2.30 


0.34 


8.81 


S23^ 


17.99 


20.85 


9.91 


17.93 


0.57 


1.20 


0.29 


8.72 


S24 


18.06 


18.59 


10.81 


22.41 


12.6 


0.44 


0.11 


8.33 


S26 


18.87 


20.06 


10.23 


17.73 


0.75 


0.50 


0.27 


8.01 


S27 


18.51 


19.46 


10.46 


21.17 


4.78 


0.92 


0.36 


8.10 



Note. — Col. (1): Target ID. Col. (2): Total extinction corrected F775W AB magnitude. 
Col. (3): Spheroid extinction corrected F775W AB magnitude. Col. (4): Logio of spheroid 
luminosity in rest frame B (solar units), not corrected for evolution. Errors are estimated to be 
0.2 d ex. Co l. (5): Spheroid effective surface brightness in rest frame B (AB magnitudes arcsec"^; 
see ^ 12.2.21 for error discussion). Col. (6): Spheroid effective radius (kpc; see ^ 12.2.21 for error 
discussion). Col. (7): Nuclear luminosity at 5100 A (10^^ erg s"-*^). Errors are estimated to be 
20%. Col. (8): Nuclear light fraction in F775W. Errors are estimated to be 20%. Col. (9): Logio 
of Mbh (solar units). Random errors are estimated to be 0.4 dex. Systematic errors are discussed 
extensively in Section |4] 

^ Spheroid size and luminosity are upper limits. 



Cosmic evolution of the Mbh — Lb relation 



TABLE 3 

Properties of the local comparison sample. 



Name 


z k 


)g Lb/Lq^b 


log Mbh /M0 


(1) 


(2) 


(3) 


(4) 


Arkl20 


0.032 


10.82 


8.18 


Mrk79 


0.022 


9.79 


7.72 


MrkllO 


0.035 


9.85 


7.40 


Mrk590 


0.026 


10.40 


7.68 


Mrk817 


0.031 


10.49 


7.69 


NGC3227 


0.004 


8.85 


7.62 


NGC4051 


0.002 


8.43 


6.28 


NGC4151 


0.003 


9.49 


7.66 


NGC5548 


0.017 


10.53 


7.83 



Note. — CoL (1): Target ID. CoL (2): Red- 
shift. CoL (3): Logio of spheroid Luminosity in 
rest frame B (solar units). Col. (4): Logio of 
Mbh (solar units). From (Peterson et al. 2004) and 
(,Bentz et al..,,2006b.) . Errors are estimated to be 0.4 
dex. 



